> 



Buchdahl-like transformations for perfect fluid spheres 

Petarpa BoonserirQ and Matt Vissei0 
School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
PO Box 600, Wellington, New Zealand 

(Dated: 2 July 2007; lOTEX-ed February 1, 2008) 



In two previous articles [Phys. Rev. D71 (2005) 124307 |gr-qc/0503007[ |, and |gr-qc/06070011 we have 
discussed several "algorithmic" techniques that permit one (in a purely mechanical way) to generate 
large classes of general relativistic static perfect fluid spheres. Working in Schwarzschild curvature 
coordinates, we used these algorithmic ideas to prove several "solution-generating theorems" of 
varying levels of complexity. In the present article we consider the situation in other coordinate 
systems: In particular, in general diagonal coordinates we shall generalize our previous theorems, in 
isotropic coordinates we shall encounter a variant of the so-called "Buchdahl transformation" , while 
in other coordinate systems (such as Gaussian polar coordinates, Synge isothermal coordinates, and 
Buchdahl coordinates) we shall find a number of more complex "Buchdahl-like transformations" and 
"solution-generating theorems" that may be used to investigate and classify the general relativistic 
static perfect fluid sphere. Finally by returning to general diagonal coordinates and making a suitable 
ansatz for the functional form of the metric components we place the Buchdahl transformation in 
its most general possible setting. 
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, I. INTRODUCTION 

Understanding perfect fluid spheres in general relativity is an important topic, recently deemed worthy of a brief 
chapter in an updated edition of the premier book summarizing and surveying "exact solutions" in general relativity [l| . 
For considerably more details, the extensive review in is an essential reference. Physically, perfect fluid spheres 
are a first approximation to realistic models for a general relativistic star. The studyof perfect fluid spheres is a 
long-standing topic with a venerable history [H, B, [p, [1, 3 , and continuing interest |l.[9l.Tl0l.[l]|. Nevertheless the 
^ ^ . topic continues to provide surprising and novel results Tl^fisl [Tsl [l^ . In particular, as derived in references 

■ and Ts*! , and as further described in references [ij, [l^ ^^ , we have developed several "algorithmic" techniques that 
I permit one to generate large classes of perfect fluid spheres from first principles in a purely mechanical way. We now 

■ generalize these algorithmic ideas, originally derived for Schwarzschild curvature coordinates, (and to a lesser extent 
^ ' isotoropic coordinates), to a number of other coordinate systems. This sometimes leads to much simpler results, and 

• ^ ■ sometimes more general results. In this current article we shall aim to: 

X ■ . 

^ ' 1. Derive a number of new and significant "transformation theorems" and "solution generating theorems" that 

5^ , allow us to map perfect fluid spheres into perfect fluid spheres. 

2. Investigate regularity conditions at the centre of the fluid sphere from the geometrical, physical, and purely 
mathematical points of view. 

3. Develop a coherent picture of how these various theorems and various coordinate systems inter-relate to one 
another. 
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II. STRATEGY 



We start the discussion by considering static spherically symmetric distributions of matter, which implies (purely 
by symmetry) that in orthonormal components the stress energy tensor takes the specific form 
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If the matter is a perfect fluid, then in addition we must have 

Pr = Pt- 



(1) 



(2) 



Invoking the Einstein equations, this then implies a differential relation on the spacetime geometry arising from 
equating the appropriate orthonormal components of the Einstein tensor 



Equivalently, one could work with the appropriate orthonormal components of the Ricci tensor 

Rqq — R'rr 



(3) 
(4) 



In terms of the metric components, this now leads to an ordinary differential equation [ODE], which constrains the 
spacetime geometry for any general relativistic static perfect fluid sphere. 

This equation constraining the spacetime geometry of static perfect fluid spheres is now analyzed in several different 
coordinate systems: general diagonal coordinates, Schwarzschild curvature coordinates, isotropic coordinates, and 
lesser-known coordinate systems such as Gaussian polar coordinates, Synge isothermal coordinates, and Buchdahl 
coordinates. In the various coordinate systems we discuss below, we exhibit the metric (in the form of the line 
element), the orthonormal components {Gff, G^g, and and Gj^) of the Einstein tensor, and the ODE resulting from 
the pressure isotropy condition Gfr = Ggg. We then classify this ODE as to its type, and whenever possible, derive 
suitable transformation theorems and solution generating theorems. 

To place the use of "unusual" coordinate systems in perspective, one might observe that Finch and Skea 0] estimate 
that about 55% of all work on fluid spheres is carried out in Schwarzschild curvature coordinates, that isotropic 
coordinates account for about 35% of related research, and that the remaining 10% is spread over multiple specialized 
coordinate systems. We take the viewpoint that the "unusual" coordinate systems are useful only insofar as they 
enable us to obtain particularly simple analytic results, and the central point of this article is to see just how much 
we can do in this regard. We also consider what conditions are required on the spacetime metric and transformation 
theorems to guarantee regularity and physicality of both the geometry and the matter distribution near the center of 
the configuaration. 



III. GENERAL DIAGONAL COORDINATES 



We begin by setting the notation. Choose coordinates to put the metric into the form: 

dr2 



A brief calculation yields 



and 



ds' = -dry dt 



RirYdn^ 



B{r) 

1 - B{R'f 2BC'R' 



i?2 



RC 



_ 1 2BC,'R' + 2CR"B + CB'R' + 2RCB + RC'B' 



_ 1 - B{R' f 2BR" - B'R' 



R 



(5) 

(6) 
(7) 

(8) 
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A. ODEs 



If this spacetime geometry is to be a perfect fluid sphere, then we niust have Gff — Ggg. This isotropy constraint 
supphes us with an ODE, which we can write in the form: 

[R{RO']B' + 2[RR"C + R^Q" - RR'C' ~ {R'fQB + 2C = 0. (9) 

In this form the ODE is a first-order hnear non-homogeneous ODE in B{r), and hence explicitly solvable. (Though it 
must be admitted that the explicit solution is sometimes rather messy.) The physical impact of this ODE ([9]) is that 
it reduces the freedom to choose the three a priori arbitrary functions in the general diagonal spacetime metric ([5]) 
to two arbitrary functions, (and we still have some remaining coordinate freedom in the r-t plane, which then allows 
us to specialize the ODE even further). 

The same ODE ([9]) can also be rearranged as: 

[2R^B]C' + [R^B' - 2BRR']C + [2B{RR" ~ [R']^) + RB'R' + 2]C = 0. (10) 

In this form the ODE is a second-order linear and homogeneous ODE in C('')- While this ODE is not explicitly 
solvable in closed form, it is certainly true that a lot is known about its generic behaviour. 
We can also view the two equivalent ODEs ([9]) and (fTO|) as an ODE for R{r): 

[2BC]RR" + [B'C - 2BC]RR' - [2BC] {R'f + [2BC' + B'C,']R^ + 2C = 0. (11) 

Viewed in this manner it is a second-order nonlinear ODE of no discernible special form — and this approach does 
not seem to lead to any useful insights. 



B. Solution generating theorems 



Two rather general solution generating theorems can be derived for the general diagonal line element. 

Theorem 1 (General diagonal 1). Suppose we adopt general diagonal coordinates and suppose that 
{C^{r)^ B{r), R{ry\ represents a perfect fluid sphere. Define 

, ( R{r)ar) V /'C'(0 {R'{r)ar)^R{rX{r)) ^ \ 

- ' [jRmw) I' y CM • (i?'MC(.)+i?WC'(r)) 'V ■ ^''^ 

Then for all A, the geometry defined by holding ^(r) fixed and setting 

ds' = -({rf dt' + + Rirfdn' (13) 

B[r) + A(r) 

is also a perfect fluid sphere. That is, the mapping 

T^eni{X):{C,B,R}^{C,B + A,R} (14) 
takes perfect fluid spheres into perfect fluid spheres 

Proof for Theorem 1. The proof is based on the techniques used in [l^ . [l3 |. No new principles are involved and we 
quickly sketch the argument. Assume that {('(?'), i3(r), i?(r)} is a solution for equation Under what conditions 
does {(l{r),B{r),R{ry\ also satisfy equation ([9])? Without loss of generality, we write 

i?(r) = B(r) + A(r). (15) 

Substitute B{r) in equation ^ 

[R{RC)']{B + A)' -f- 2[RR"C + - RR'C' ~ {R' f(\{B + A) -f 2C = 0. (16) 

Rearranging we see 

[R{RC)']B' + 2[RR"C + R^C - RR'C - {R'fQB + 2C (17) 
+ [R{RCY]A' + 2[RR"C + R^C' - RR'C ~ {R'fC^ = 0. 
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But we know that the first hne in equation (jT?!) is zero, because {R, C, B} corresponds by hypothesis to a perfect 
fluid sphere. Therefore 

[i?(i?C)'] A' + 2[RR"C + R^C - RR'C - {R'?Q^ = 0, (18) 

which is an ordinary homogeneous first order differential equation in A. A straightforward calculation, including an 
integration by parts, now leads to 

as required. □ 

Theorem 2 (General diagonal 2). Suppose we adopt general diagonal coordinates, and suppose that 
{Q{r), B{r), R{ry\ represents a perfect fluid sphere. Define 

Z{r) =a + e f dr. (20) 

J ./mar? 

Then for all a and e, the geometry defined by holding Bir) and R{r) fixed and setting 

ds2 = ^(^{r? Z{rf At^ + — - + R^rfAV? (21) 

B{r) 

is also a perfect fluid sphere. That is, the mapping 

Tgen2{<y, e) : {C, B, R} ^ {C Z{C, B, R), B, R} (22) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof. The proof is based on the technique of "reduction in order" as used in [T^ . [T3 | . No new principles are involved 
and we quickly sketch the argument. Assuming that {C,{r), B{r), R{r)} solves equation (jlOp . write 

C(r) ^ C(r) Z{r) . (23) 

and demand that {C,{r) Z{r), B{r), R{r)} also solves equation (ITUl) . Then 

[2R^B] {(Z)" + [R^B' -2BRR'\{CZ)' + [2B{RR" - [R'f)+RB'R' + 2]{CZ) = 0. (24) 
We expand the above equation to 

[2i?2 B] iC" Z + 2C Z' + C Z") + B' ~2BR R'] {(' Z + ( Z') 

+ [2B{RR" ~ [Rf) + RB'R' + 2]{CZ) = 0, (25) 

and then re-group to obtain 

{[2i?2 B] C" + [i?^ B' -2BRR']C + [2B{RR" ~ [Rf)+RB'R' + 2]C} Z 

[2R^ B] {2C Z' + ( Z") + B' -2BR i?']C Z' = 0. (26) 

This is a linear homogeneous 2nd order ODE for Z . But under the current hypotheses the entire first line simplifies 
to zero — so the ODE now simplifies to 

B' C + 4 R^ B C - 2B RR' C] Z' + {2R^ B () Z" = 0. (27) 

This is a first-order linear ODE in the dependent quantity Z'(r). Rearrange the above equation into 

Z" [R^B'C + AR^BC -2BRR'C] 



Z' (2i?2 B C) 

Simplifying 



(28) 



z" 1 B' c R' 

-W--2B-'-C^R- (2^) 
Re-write Z" jZ' — (lln(Z')/dr, and integrate twice over both sides of equation (PH)). to obtain 

^ dr, (30) 

depending on the old solution {C, B, i?}, and two arbitrary integration constants a and e. 

□ 
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C. Regularity conditions 



Now consider the result of placing various regularity conditions on the geometry of the perfect fluid sphere. (See 
particularly the discussion by Delgaty and Lake regarding the importance of regularity conditions At a minimum 
we want the spacetime geometry to be well-defined and smooth everywhere. This requires: 

• (^{r) > 0, so that a clock held at any fixed but arbitrary "place" (r, 9, (f) is well behaved, and at least continues 
to "tick" (even if it is fast or slow). 

• Now the "centre" of the spacetime is at some a priori arbitrary coordinate tq such that the very specific condition 
R{ro) = is satisfied; by a simple coordinate change we can without loss of generality set r-g ^ so that at the 
centre of the spacetime R{0) — 0. 

• So purely for geometrical reasons we should demand: 

- i?(0) = 0; 

- C(0) > 0. 

• Now consider the spatial three-geometry 

dr^ 

+ Rir)'dn\ (31) 

B{r) 

and consider a small "ball" centered on the origin with coordinate radius r. The surface area of such a ball is 

S = A7rR{rf, (32) 

while its proper radius is 

dr 

(33) 



VW) 

If the geometry at the centre is to be smooth (that is, locally Euclidean) then these two measures of radius must 
asymptotically agree (for small balls): 

t = R + 0{R^). (34) 

Thus in terms of the original coordinate r 

dR 



di 
dr 



r=0 VW) 



(35) 

r=0 



implying 

B{0) [i?'(0)]2 = 1, and B{0) ^ 0, that is B{0) = 1/[R'{0)]'^. 

• Now if the spacetime is regular at the centre we must avoid "kinks" in the metric components. Specifically, 
in terms of proper distance from the centre we should have (,{£) — C(0) -I- 0{£'^). But in view of the already 
established relation between £ and r this imphes ({r) = C(0) + 0{r'^) and we must have 

C'(0) = 0. (36) 

• This does not complete the list of regularity conditions. We can also look at the orthonormal components of 
the Einstein tensor. Consider 

From the conditions already established, applying the I'Hospital rule to the last term leads to a finite result 
2i3(0)C"(0)/C(0). For the remaining term, noting that the denominator (i?^) has a double zero at the centre. 
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we deduce that the derivative of the numerator must hkewisc be zero, in order that there be any chance for the 
I'Hospital rule to dehver a finite answer. This now yields the additional constraint that at the centre 

[B{R'f]' -y 0, (38) 

that is 

B'(0)J?'(0) + 2S(0)i?"(0) = 0. (39) 
In view of our earlier result for B(0) this can be written as 

2i?"(0) / 1 



2i?"(0) _ / 



[i?'(0)]3 \[R'{r)f 
To check consistency, consider 



(40) 



r=0 



a„ + c,,'"''^''"^- 1" (41) 



and note that this is finite at the centre provided C' ^ and {B[R'Y\' 0. 

• As a final consistency check consider 

_BCR' , 2C"B + CB' , 2R"B + B'R' 
^ee - - + + ^ , (42) 

and again note that this is finite at the centre provided (' ^ and [JB(i?')^]' — + 0. 

• So now let us collect all our regularity conditions: 

RiO) = 0, C(0)>0, B(0) = -^>0, C'(0) = and B'(0) = -^|^. (43) 

• Finally inserting all these conditions into the r — > limit we quickly find that at the centre of the fluid sphere 
the Einstein tensor has the limit: 

r I - r \ -r \ - ^ + 2^(0) ^ R"'{0) 3[R"m' 
^ff\o - ^ee\o-^U\o- 2 + C(0) [i?'(0)]2 + [i?'(0)]3 [i?'(0)]4 ' 



3 [p"rm i ^^'"^^^ m"m' ' 

2 \ ^ [R'{0W [-R'(o)]^ i ' 



Gu\ o = --^\ B"(0) + - \ , (45) 



and therefore 



r \ ^-^r I 6C^'(0) 
^"'° + ^^^^l° = C(0)[i?'(0)]2- 

• There is another largely independent way of getting these regularity conditions, either by hand or with the aid 
of a symbolic manipulation program such as Maple. Let us start by assuming that near the origin all metric 
functions have Taylor series expansions 

C(r) =Co + Cir+^ + ^+ Oir% (47) 
B{r)=Bo + B,r+^ + ^ + 0{r% (48) 
Rir) =R,r+^ + ^ + Oir% (49) 
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Now use this ansatz to calculate G^f, Gff^ and Ggg as Laurent series in r. A symbolic manipulation program 
such as Maple will quickly tell you 



G,i, = ^ + ^ + C,i, + 0{r). (50) 



Regularity now has the obvious meaning that you want the coefficients A- and B^^^ of the f /r^ and f /r pieces 
to vanish — which implies constraints on the Taylor series coefficients {Q, Bi, Ri}. A brief computation should 
convince one that the pole residues A.g and -B-g vanish provided: 

i?o = 0, Co > 0, Bo = l/Rl > 0, Ci = 0, and Bi = -2i?2/(i?i)'; (51) 

that is 

i?(0) = 0, C(0)>0, B(0) = l/[i?'(0)2] > 0, C'(0) = 0, and B'(0) = -2i?"(0)/[i?'(0)]3. 

(52) 

These are the same regularity conditions as we previously found by hand. Imposing these conditions, symbolic 
manipulation software such as Maple can then easily evaluate the limits as you move to the centre of the star. 

The situation for general diagonal coordinates is the most complicated case considered in this article, other cases 
discussed below are simpler. 

D. Implications 

Note that the regularity conditions do feed back into theorem General diagonal 2. In that theorem we were 
interested in the quantity 

Z{r) =a + e f J!^^ dr. (53) 

J ./mar? 

In view of the regularity conditions we now know that for a physically reasonable fluid sphere the integrand is finite 
all the way to the origin so that we can without loss of generality set 

Z{r)=a + e( dr. (54) 

But to ensure regularity we want both ({0} and the new C(0) — ^ C(0) Z{0) to be finite and positive, hence Z{Q) > 0, 
and therefore a > 0. This reproduces in a purely mathematical way one of the results of reference ^13] , where we used 
the TOV and boundedness of the central pressure and density to impose a similar constraint. 



E. Summary 



Let us now assess what we have done so far: For arbitrary spherically symmetric static spacetimes in general 
diagonal coordinates we have investigated the perfect fluid constraint, derived two "solution generating theorems", 
and applied regularity conditions to the spacetime. The analysis has so far been rather general, and we still have the 
freedom to make one further specific coordinate choice to simplify the line element, Einstein tensor, solution generating 
theorem, and regularity conditions. Making the extra coordinate choice will sometimes simplify the mathematical 
structure sufficiently to allow more information to be extracted. 



IV. SCHWARZSCHILD CURVATURE COORDINATES 



The Schwarzchild curvature coordinate system is by far the most popular coordinate system used in the study of 
perfect fluid spheres, and this coordinate choice (corresponding to R{r) r) accounts for approximately 55% of 
the research reported on perfect fluid spheres H. The Schwarzschild coordinate system is the basis of the various 
theorems derived and discussed in references [l^Tlsl. [l3. [isl [l6t . The line element of a spherically symmetric spacetime 
in Schwarzschild curvature coordinates is 

ds2 = _((^)2 + ^ + r^an\ (55) 
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To begin, we calculate 

Gff — — ^ — — ^ — I (56) 
^§§-'^^-2 ^ ■ ^ ' 

and 

A. ODEs 

Pressure isotrropy leads to the ODE 

[r{rCy]B' + [2r^C" ' 2«)']S + 2C = 0. (59) 

This is first-order linear non-homogeneous in B(r). Solving for B{r) in terms of C(r) is the basis of the analysis in 
reference 11 1, and is also the basis for Theorem 1 in reference |lj|. (See also [ij and the reports in [l^, [l^. After 



rephrasing in terms of the TOV equation this is also related to Theorem P2 in reference [13|). If we re-group in terms 
of C(r) we find 

2r2c" + {r^B' - 2rB)C + [rB' - 2B + 2)C = 0, (60) 

which is a linear homogeneous second-order ODE. This is the basis of Theorem 2 in [l^, and after being recast as a 
Riccati equation is also the basis of Theorem PI in reference [isj . 

B. Regularity conditions 

Now consider the regularity conditions we demand to keep the spacetime regular at the centre. First we have: 

• The "centre" of the spacetime is by construction at r = 0. 

• We must again demand ((r) > 0, so that a clock held at fixed (r, 9, (f) is well behaved and at least continues to 
"tick" (even if it is fast or slow). 

• In particular ^(0) > 0. 

• We now consider the spatial geometry 

and slightly modify the discussion (previously developed in general diagonal coordinates) for a small "ball" 
centered on the origin with coordinate radius r. The surface area of such a ball is 

5 = 47rr2, (62) 

while its proper radius is 

dr 



, (63) 

/o VW) 

If the geometry at the centre is to be smooth (locally Euclidean) then these two measures of radius must 
asymptotically agree (for small balls): 

l^r + Oir"^). (64) 
But in terms of the original coordinate r this implies 



5(0)^^0, and ^^(0) = 1, that is B(0) = 1. (65) 
Note we now have the much tighter constraint -B(O) = 1. 
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• The same argument used for general diagonal coordinates still implies 

C'(0) = 0. (66) 



• By considering 



we deduce that at the centre 



That is 



To check consistency note 



G.. = ^ + ^ (67) 

B' 0. (68) 
B'{0) = 0. (69) 



G,, + G,, = ^^-^^^, (70) 

which is finite at the centre provided B' — > and — > 0. 

• For the final consistency check note 

^ee -'-4>4>-^ + ^ + 7^' (71) 

which is finite at the centre provided ^ and C' ^ 0. 

• So now let's collect all our regularity conditions: 

C(0) > 0, B{0) = 1, C'(0) = and B'{0) = 0. (72) 

• Inserting all these conditions into the limit we quickly find 

r \ -r \ -r \ _ 2C"(o) , b"(o) 

Wtio- ^ — , (74) 



and 



(jfflo + 3(jrf|0 - ^^Q^ • (75) 

• The regularity regularity conditions derived above can again also be extracted in another manner — by using a 
Taylor series expansion for the metric components to derive a Laurent series for the Einstein tensor — using (for 
example) Maple, this will still give the same regularity conditions as we derived by hand. This is a specialization 
of the method developed in the previous general diagonal coordinates. A little work with Maple should convince 
you that the pole pieces vanish provided: 

C(0) > 0, B{0) = 1, C'(0) = and B\0) = 0. (76) 

The ultimate consistency check on this formalism now results from considering general diagonal coordinates and 
simply setting R{r) r. That the results agree under this substitution gives us confidence that formalism is 
internally consistent. 
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V. ISOTROPIC COORDINATES 

Isotropic coordinates are again commonly used coordinates when investigating perfect fluid spheres. This is the 
second most popular coordinate system, accounting for about 35% of published research in this field [sj. The spacetime 
metric for an arbitrary static spherically symmetric spacetime in isotropic coordinates is conveniently given by 

d.^ = -Cirf dt' + ^^-^^L_{d.2 + rW}. (77) 

We calculate 

Gff = -2BB' CVr + (B'fC^ - C'^B^, (78) 
G§e = = -BBXVr + (B'f ^ ~ BB''^ + B^C'\ (79) 

and 

= 2B'^(C" + ^B'^C.C/r ~ 3B^C^ - 2BB'CC + 2BB"C^ ~ 3B'^C^ + iBB'C^/r. (80) 

A. ODEs 

The pressure isotropy condition leads to the very simple looking ODE 0: 

'(^'Y B"-B'/r 



C / 2B 



(81) 



There are several ways of improving this. For instance, if we write C(^) = exp(J (7(r)dr) then we have an algebraic 
equation for ^(r) [9,]: 



/ B" - B'/r 
g(0=±V 2B ■ 

Conversely, the isotropy condition can be written in terms of B{r) as: 

B" - B'/r - 2g^B = 0. (83) 
There is also an improvement obtained by writing B(r) — exp(2 J h{r)dr) so that 

g{r)^ ^2h{r)^ + h'{r)-h{r)/r (84) 
which is the basis of the analysis in reference [§] . 

B. Solution generating theorems 

Let us now develop two transformation theorems appropriate to isotropic coordinates. 

Theorem 3 (Isotropic 1 — Buchdahl transformation). In Isotropic coordinates, if {({r), B{r)} describes a 
perfect fluid then so does {C,{r)^^ , B{r)}. This is the well-known Buchdahl transformation in disguise. That is, if 

ds' = -dry dt' + ^^-^J^{dr2 + r'dn'} (85) 



represents a perfect fluid sphere, then the geometry defined by 

d.^^-^dt^ + |!i^{dr^+rW} (86) 



is also a perfect fluid sphere. Alternatively, the mapping 

Tisoi:{C,B}^ {C\B} (87) 
takes perfect fluid spheres into perfect fluid spheres, and furthermore is a "square root of unity" in the sense that: 

Tiso 1 ° Tiso 1=1- 
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Proof. By inspection. Vide equation ([81]) . □ 
Theorem 4 (Isotropic 2). Let {({r), B{r)} describe a perfect fluid sphere. Define 

rdr 

mr) 



Z={<^ + ^1^\- (89) 



Then for all a and e, the geometry defined by holding C(r) fixed and setting 

= ~ar? + ^^-^^-L__{,ir^ + r'an'y (90) 

is also a perfect fluid sphere. That is, the mapping 

Tiso2i^,e):{C,B}^{C,B Z{B)} (91) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof. The proof is based on the technique of "reduction in order" , and is a simple variant on the discussion in . 
Assuming that {({r), B{r)} solves equation (|83l) . write 

B{r) ^ B{r) Z{r) . (92) 

and demand that {(^{r), B{r)Z{ry\ also solves equations (fSTj) and ((83|) . We find 

{B Z)" -{B Z)' /r-2g^{B Z) = Q. (93) 

Re- grouping 

{B" - B' /r - 2g'^B] Z + {2B' - B/r)Z' + BZ" = 0. (94) 

This is a linear homogeneous second-order ODE for Z which now simplifies [in view of ((83|) ] to 

(2B' -B/r)Z' + BZ" ^0, (95) 

which is an ordinary homogeneous second-order differential equation, depending only on Z' and Z" . (So it can be 
viewed as a first-order homogeneous order differential equation in Z', which is solvable.) Separating the unknown 
variable to one side, 

Z' Br 

Re-write Z" jZ' = dln(Z')/dr, and integrate twice over both sides of equation (|96p. to obtain 

depending on the old solution {((r), B{r)}, and two arbitrary integration constants a and e. □ 
Theorem 5 (Isotropic 3). The transformations Tiso i and Tiso 2 commute. 

Proof. By inspection. □ 

Note that the fact that these two transformation theorems commute is specific to isotropic coordinates. Such 
behaviour certainly does not occur in Schwarzschild coordinates or in general diagonal coordinates. 
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C. Regularity conditions 



We now investigate the regularity conditions for isotropic coordinates. Some details are a little tricky, and have led 
to confusion in some of the previous published literature. 



• The "centre" of the spacetimc is at some coordinate ro; by a simple coordinate change we can set ro 

that at the centre of the spacetimc is at r = 0. 



so 



• We still require C(r) > 0, so that a clock held at fixed (r, 9, (p) is well behaved and at least continues to "tick" 
(even if it is fast or slow). In particular 



C(o) > 0. 



Now consider the spatial geometry 



dr^ 



+ ■ 



(98) 



(99) 



C(r)2B(r)2 C{rfB{r)^ C{rf B{rf 

and consider a small "ball" centered on the origin with coordinate radius r. The surface area of such a ball is 

^2 

477 (100) 



while its proper radius is 



C(r)2 S(r)2 ' 



dr 

ar)B(ry 



(101) 



If the geometry at the centre is to be smooth then these two measures of radius must asymptotically agree (for 
small balls): 



dr 



Q{r)B{r) C(r-)S(r) 



implying 



that is 



ar)B{r) \ar)B{r) 



+ 0{r% 



0(r), 



C(r)B(r) C{r)B{r] 



+ 0{r). 



(102) 



(103) 



(104) 



This implies no constraint on (^(r) or B{r), beyond the obvious constraint that these quantities remain finite 

and nonzero. Note in particular that in isotropic coordinates it is not true that we need to enforce -B(O) = 1. 
(If B(0) 1 we could think of arbitrarily rescaling it by a constant factor to force it to be unity, this would at 
worst multiply all curvature invariants by a constant, and so a nonsingular geometry remains nonsingular under 
such rescaling. This point is commonly misunderstood in the literature.) 

On the other hand, we could also find all our regularity conditions for isotropic coordinates by deducing them 
from the analysis developed for general diagonal coordinates. Identifying coefficients in the two metrics we see 



c 



gen i 



iCB) 



2 

iso 



B 



gen-i 



and 



Now consider 



{QB), 



R 



gen- 



[B{R'f 



(cs)io 



r_C 
C 



rB' 



(105) 
(106) 

(107) 



So at the centre [B{R'Y]gp.a 1 provides no extra information about Bjso or (beyond the fact that these 
quantities are finite and nonzero). 
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• The same argument used for general diagonal coordinates and Schwarzschild coordinates still implies 

C'(0) = 0. (108) 

• This does not complete the list of regularity conditions. Consider 
Finiteness of this quantity now implies 

B'{0) = (110) 

• On the other hand, finiteness of 

Gee = = + B'^C - BB"e + BX', (m) 

implies no additional constraint. 

• Similarly, finiteness of 

Gff = 2B\C + - SB'C" - 25B'CC' + 2BB"C'' - SB'^C' + ^ , (112) 

implies no additional constraint. 

• Summarizing, our regularity conditions are 

B{0) > 0; B'{0) = 0; C(0) > 0; C'(0) = 0. (113) 

• Now inserting the above conditions into the central limit we quickly find 

''eel 



Grf\ = Ggs\o = G, ,| = -2B(0)S"(0)C(0)^ (114) 



Gfil = 6B(0)B"(0)C(0)2 + 6B{0fmC"{0), (115) 

and 

GhIo + 3G,Ho = 6B(0)2C(0)C"(0). 

• On the other hand, the regularity conditions can also be easily obtained by using Taylor series for the met- 
ric components to compute Laurent series for the orthonormal components of the Einstein tensor. A brief 
computation with Maple should convince you that the pole pieces vanish provided: 

C(0) > 0, B{0) > 0, C'(0) = and B'{0) = 0. 
D. Implications 

Note that the regularity conditions do feed back into theorem Isotropic 2. In that theorem we were interested in 
the quantity 

^W=- + e/;^. (116) 

In view of the regularity conditions we now know that for a physically reasonable fluid sphere the integrand is flnite 
all the way to the origin so that we can without loss of generality set 

But to ensure regularity we want both B{Q) and the new B{0) — > B{Q) Z{0) to be finite and positive, hence Z{0) > 0, 
and therefore a > 0. 
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VI. GAUSSIAN POLAR COORDINATES 

Having now dealt with the most common coordinate systems, viz. Schwarzschild curvature coordinates and isotropic 
coordinates, we wiU spend some time studying other more unusual coordinate systems. Although Gaussian polar 
coodinates are typically less convenient than Schwarzschild or isotropic coordinates, we have however found them to 
be occasionally useful in the theory of symmetric static perfect fluid spheres. Consider, for instance, the metric in 
Gaussian polar (proper radius) coordinates: 

ds^ = -CW^ di^ + dr^ + R{r)^dn^. (118) 

The following expressions can easily be calculated 

a,,,a,,,<J^^^^, (120, 

and 

Cu ^ -^"-^f""' . (121) 

The isotropy condition supplies a second-order linear homogeneous ODE for C{r): 

In contrast, when rearranged and viewed as an ODE for R{r), 

QRR" - C'RR' - CiR'f + C'R^ + C = 0, (123) 
this has no particularly desirable properties. 

A. Solution generating theorem 

Theorem 6 (Gaussian polar). Suppose we are in Gaussian polar coordinates, and that {C(r),i?(r)} represents a 
perfect fluid sphere. Define 

Then for all a and e, the geometry defined by holding R{r) fixed and setting 

ds^ = ^C{r)^A{rf dt^ + dr2 + R{r)^dn^ (125) 
is also a perfect fluid sphere. That is, the mapping 

TcsussisnicT, e) : {C, R} ^ {C A(C, R), R} (126) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof. The proof is again based on the technique of "reduction in order" . No new principles are involved and the 
argument is quickly sketched. Assuming that {({r), R{r)} solves equation (|122p . write 

C(r) ^ CM A(r), (127) 

and demand that {C{r) A(r), R{r)} also solves equation ()122|) . Then 

(C A)" - (C A)'^ + (C A) |W^T±^| ^ 0^ (128) 
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which simplifies to 



2C'-C^) A' + CA" = 0. (129) 



This is a Ist-ordcT homogeneous order differential equation in A', which is solvable. Separating the unknown variable, 
and integrating twice one obtains 

AW = f + ./^}- ("0) 

depending on the old solution {C{r), R{r)}, and two arbitrary integration constants a and e. □ 

B. Regularity conditions 

In the following section, wc briefly describe the regularity conditions that must hold at the centre of the fluid sphere 

• Adapting the arguments that by now should be completely standard we have 

i?(0) = 0; C(0)>0; C'(0) = 0. (131) 

• Now consider the spatial geometry 

dr^ + i?(r)2 dO^ (132) 
and consider a small "ball" centered on the origin with coordinate radius r. The surface area of such a ball is 

S" = 47r i?(r)2, (133) 

while its proper radius is 

1 dr = r. (134) 





If the geometry at the centre is to be smooth, these two measures of radius must asymptotically agree (for small 
balls): 

r = R + 0{R'^). (135) 
But in terms of the original coordinate r this implies 

i?'(0) = 1. (136) 

Now consider 

_ R'^ 1 2CR'/C 

The last term is finite by the I'Hospital rule and the previously derived condition 0. Applying the I'Hospital 
rule twice the first term will be finite at the centre iff 

[{RTy - 0, (138) 

that is, iff 

R"{0) = 0. (139) 

For a consistency check, note that 

G,, + Gu=^-^^^^, (140) 

which implies no new constraint. 
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For a consistency check, note that 



which impHes no new constraint. 

• Now let us collect all our regularity conditions: 

R{0) = 0, R'{0) = 1, i?"(0) = 0, C(0) > 0, and C'(0) = 0. (142) 

• Inserting all these conditions into the limit we quickly deduce 

Grf\o = G,-^| = G^^l = R"'{0) + (143) 

G,-,-|o = -3i?"'(0), (144) 

and 

Gttio + JGfflo - ^^^^ ■ (145) 

• These regularity conditions can again be verified by using Taylor series in the metric components to calculate 
Laurent series for the orthonormal components of the Einstein tensor. Using Maple, the pole pieces vanish 
provided: 

C(0) > 0, C'(0) = 0, i?(0) = 0, R'{0) = 1, and R"{0) = 0. (146) 

C. Implications 

Note that the regularity conditions do feed back into theorem Gaussian polar. In that theorem we were interested 
in the quantity 

In view of the regularity conditions wc now know that for a physically reasonable fluid sphere the integrand is finite 
all the way to the origin so that we can without loss of generality set 



Z{r) =a + e [ 
Jo 



But to ensure regularity we want both ^(0) and the new ^(0) ({0) Z{0) to be finite and positive, hence Z{0) > 0, 
and therefore a > 0. 



VII. SYNGE ISOTHERMAL COORDINATES 

Consider the metric in Synge isothermal (tortoise) coordinates: 

ds2 = -C(r)-2 {df - dr^} + {C(r)-2 R{rfdn^}. (149) 

A brief computation yields 

3R^ icr - mR c + {Rre - e 

Lrff — ^2 ' \^'^^) 



3R{Cf-2CR'C-2RCC" + eR" 
'^ee ^ ' \^^^) 
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and 

Git = • (152) 

Pressure isotropy supplies us with an second-order homogeneous ODE for C(^): 

,-._C.g_; p-W';""-- Uo, (153, 



R \ 2i?2 

This is very similar to the isotropy equation in Gaussian polar coordinates. On the other hand, rearranging to obtain 
an ODE for R{r) yields a second order nonlinear difFerntial equation of no particularly useful form. 

A. Solution generating theorem 

Our first key result can be phrased as a simple theorem: 
Theorem 7 (Synge). Suppose {C(r),i?(r)} represents a perfect fluid sphere. Define 

Then for all a and e, the geometry defined by holding R{r) fixed and setting 

is also a perfect fluid sphere. That is, the mapping 

Tsyngei<y, c) : {C, R} ^ {C ^(C, R) , R} (156) 

takes perfect fluid spheres into perfect fluid spheres. 

Proof. The proof is again based on the technique of "reduction in order" . No new principles are involved and the 
argument is very quickly sketched. Assuming that {C(r),i?(r)} solves equation (jl53[) . write 

C(r) CM Mr) ■ (157) 
and demand that {C{r) A{r), R{r)} also solves equation (|153p . We eventually find 

2C'-C ^V' + C^" = 0, (158) 



which is solvable. Separating the unknown variable to one side, and integrating twice we obtain 



R{r) dr 



A{r)^{a + e j ^^1^ \ . (159) 



depending on the old solution {((r), R{r)}, and two arbitrary integration constants cr and e. □ 

B. Regularity conditions 

• As usual 

i?(0) = 0; C(0) > 0; C'(O) = 0. (160) 

• Now consider the spatial geometry 

dr2 , R{r)^ ^ dr^ + Rjrf dn^ ^^^^^ 
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and consider a small "ball" centered on the origin with coordinate radius r. The surface area of such a ball is 

^ = 4-^, (162) 

while its proper radius is 

1 



dr. (163) 



c(r-: 

If the geometry at the centre is to be smooth then these two measures of radius must asymptotically agree (for 
small balls): 

e=^ + 0[iR/Cn (164) 

But in terms of the original coordinate r this implies 

R'{0) = 1. (165) 
Note that there is no additional constraint on ^(0) apart from the fact that it is nonzero and finite. 
Consider 

G,, = SR'^ - ^ + ^-^^ (166) 

Regularity at the centre requires 

mTr 0, (167) 

which implies 

i?"(0) = 0. (168) 

The combination 

Grr + Gu = 2a"-^^, (169) 

does not lead to any new constraint. 
The component 

Gee = % = 3(C')^ - 2CC" + ^^^'^^ (170) 

does not lead to any new constraint. 
Collecting all our regularity conditions: 

C(0) > 0, C'(0) = 0, i?(0) = 0, i?'(0) = 1, and i?"(0) = 0. (171) 

Finally inserting all these conditions into the limit we quickly find 

Grflo = Ggglo = G^^lo = -4C(0)C"(0) + C(0)^i?"'(0), (172) 

GhIo = 6C(0)C"(0) -3C{0fR"'{0), (173) 

and 

G,-,-|o + 3G,^|o = -6C(0)C"(0). (174) 

Again working with Taylor series to derive Laurent series, the pole pieces vanish and the centre of the sphere is 
regular provided: 

C(0) > 0, C'(0) = 0, i?(0) = 0, i?'(0) = 1, and i?"(0) = 0. (175) 
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C. Implications 

Note that the regularity conditions do feed back into theorem Synge. In that theorem we were interested in the 
quantity 

In view of the regularity conditions we now know that for a physically reasonable fluid sphere the integrand is finite 
all the way to the origin so that we can without loss of generality set 

r R{r)dr 

But to ensure regularity we want both ({0) and the new ^(0) ({0) Z{0) to be finite and positive, hence Z{0) > 0, 
and therefore ct > 0. 

VIII. BUCHDAHL COORDINATES 

Without loss of generality we can put the metric in "Buchdahl coordinates" and choose the coefficients to be 

ds2 = -(^{rf dt^ + C(r)-2 {dr^ + R{rf dn^} . (178) 

This coordinate system is a sort of cross between Synge isothermal (tortoise) coordinates and Gaussian polar (proper 
radius) coordinates. We calculate 

G,, = -iCr- ^" (179) 
Gee = Gu = +iCr + ^-^, (180) 

and 

G,, = -3,cV + 2C("-M5!^^ClE_(im. ,181) 
Imposing pressure isotropy supplies us with a particularly simple first-order homogeneous ODE for C{r): 



C J 2i?2 ■ 

This is very similar to the equation we obtained in isotropic coordinates. (Rearranging this into an ODE for R{r) 
yields a second-order nonlinear differential equation which docs not seem to be particularly useful.) 

Theorem 8 (Buchdahl). If {({r), R{r)} describes a perfect fluid then so does {({r)~^ , R{r)} . This is the Buchdahl 
transformation in yet another disguise. The geometry defined by holding R{r) fixed and setting 

ds^ = -C(r)-' dt^ + C(r)2 {dr' + R{rfdn''} (183) 

is also a perfect fluid sphere. That is, the mapping 

TBuchdahi:{C,R}^ {C\R} (184) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof. By inspection. (Note strong similarities to the discussion for isotropic coordinates.) □ 
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A. Regularity conditions 

As usual 

i?(0) = 0; C(0) > 0; C'(0) = 0. (185) 

Now consider the spatial geometry 

C(r)-2 {dr^ +i?(r)2dfi2} (186) 
and consider a small "ball" centered on the origin with coordinate radius r. The surface area of such a ball is 

S = 4TTC{r)-^ R{rf, (187) 

while its proper radius is 



£■■ 



[ C(r)"Mr. (188) 
Jo 



If the geometry at the centre is to be smooth then these two measures of radius must asymptotically agree (for 
small balls): 

l = R/C + 0[{R/Cf]. (189) 
But in terms of the original coordinate r this implies 

R'{Q) = 1 (190) 

Consider 

G.. = -(CT-^^i^^P^ (191) 

this now yields the additional constraint 

[{R!f]' ^ 0, (192) 

which implies 

R"{Q) = 0. (193) 

The combination 

G,, + = -4{Cr + 2CC" - 2C(Ci?"-2C'i?0 ^ (194) 

does not lead to any new constraint. 
Similarly 

Gee = iCr + ^-^, (195) 

does not lead to any new constraint. 

So now let's collect all our regularity conditions: 

R{0) = 0, R'{0) = 1, R"{0) = 0, C(0) > 0, and C'(0) = 0. (196) 

Finally inserting all these conditions into the limit we quickly find 

Gfr\o = G^llo = G^^lo = 2C(0)2 i?"'(0), (197) 

Gii\ = 6C(0)C"(0) - 6C(0)2 R"'{0), (198) 

and 

Gjj|o + 3G,,|o = 6C(0)C"(0). (199) 
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• An analysis based on Taylor series and Laurent series again quickly leads to : 

i?(0) = 0, R'{0) = 1, R"{0) = C(0) > 0, and C'(0) = 0. (200) 

We have now seen the Buchdahl transformation ({r) — > l/Cir) arise in two separate settings — in both isotropic 
coordinates and in the Buchdahl coordinates we have just investigated. This raises the question as to whether we 
can expect to encounter Buchdahl like transformations in more general settings, perhaps by suitably choosing the 
functional form of the metric. This will be our last topic to be investigated in this article. 



IX. GENERALIZED BUCHDAHL ANSATZ 



We now return to general diagonal coordinates and, based on our insight from dealing with isotropic coordinates 
and Buchdahl coordinates, make a specific ansatz for the functional form of the metric components. Without loss of 
generality we choose the metric coefficients to be 

ds2 = -C(r)2 dt^ + C(r)-2 \^^+ RirfdQ''^ . (201) 



We calculate 



G,. = -EiCr - (202) 



and 



Imposing pressure isotropy supplies us with a first-order homogeneous ODE for C(^): 

'C'V [2-2E{R')'^ + 2ERR" + RE'R'] 



C J 4:ER'^ 



(205) 



This is very similar to the equation we obtained in isotropic coordinates and Buchdahl coordinates, but now in general 
diagonal coordinates — the key point is that we have carefully chosen the functional form of the metric components. 
Rearranging this into an ODE for E{r) yields a first-order linear differential equation 

[C'^RR']E' + [4i?2((')2 + 2RC^R" - 2(^{R'f]E + 2(^ ^ 0. (206) 

(In contrast, rearranging this into an ODE for R{r) yields a second-order nonlinear differential equation which does 
not seem to be particularly useful.) 

Theorem 9 (Generalized Buchdahl). // {(^(r), £'(r), i?(r)} describes a perfect fluid then so does 

{C,(r)~^ , E(r), R{r)}. This is the Buchdahl transformation in yet another disguise (now in its most general setting). 
The geometry defined by holding E(r) and R{r) fixed and transforming the given perfect fluid sphere 

„2 



into the new geometry 



ds^ = -Cirf dt^ + C(r)-2 Rirfdn''^ • (207) 



ds^ = -C(r)-2 dt^ + C(r)2 \^-^+ RirfdO:' } . (208) 



yields a spacetime that is also a perfect fluid sphere. That is, the mapping 

"^Generalized Buchdahl • 

{C,E,R}^ {C\E,R} (209) 

takes perfect fluid spheres into perfect fluid spheres. 
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Proof. By inspection. (Note very strong similarities to the discussion for isotropic coordinates and Buchdahl coordi- 
nates. Note that because we are now working in general diagonal coordinates, where the metric has three independent 
components, this is the most general result we can hope for for a static spherically symmetric spacetime.) □ 

The regularity conditions in this case can be extracted from a minor variant of the previous analysis, and reduce 
to: 

i?(0) = 0, i;(0)[i?'(0)]2 = 1, £;'(0) = -2i?"(0)/[i?'(0)]3 = 0, C(0)>0, and C'(0) = 0. (210) 
Finally inserting all these conditions into the central limit we quickly find 

^1 ^1 ^1 C(o)^ R"'io) K{orR"io? , C(o)^i^"(o) 

Cxff I - OggI - G^^l - -^7(^ RUfij^ 2 ' ^ ^ 



_ -3m'R"'iO) , 9C(0)^fi"(0)^ 3C(0)^i?"(0) , 6C(0)C"(0) 

i?'(0)3 ^ i?'(0)4 2 i?'(0)2 ' ^^'^^ 

and 

Gjjl + 3 Gff 1 = ^jJ^Qy ^ ■ 

In addition, we could use the fact that the isotropy condition, when viewed as a differential equation in E{r), is 
first-order linear to develop yet another variant on the theorem General diagonal 1. As no new significant insight 
is gained we suppress the details. 



X. DISCUSSION 



In this article, we have developed several new and significant transformation theorems that map perfect fluid spheres 
to perfect fluid spheres using both "usual" and "unusual" coordinate systems — such as Schwarzschild (curvature), 
isotropic, Gaussian polar (proper radius), Synge isothermal (tortoise), and Buchdahl coordinates. In each case we 
developed at least one such transformation theorem, while in several cases we have been able to develop multiple 
transformation theorems. 

We have also investigated regularity conditions at the centre of the fluid sphere in all these coordinate systems. 
The meaning of regularity is that (at a minimum) we want the spacetime geometry to be well-defined and smooth 
everywhere. We can calculate the regularity conditions by hand, but note that there is another way to more me- 
chanically derive these regularity conditions by using Taylor and Laurent series, with calculations aided by symbolic 
manipulation programs such as Maple. We have found it convenient to use this observation to provide a crosscheck 
on our explicit computations. 

In summar y, w e have now extended the algorithmic technique originally introduced in reference 9] and extended 
in references |1Q| and [ll|, [l^, [H, [13, [iM EBl to many other coordinate systems and many other functional forms for 
the metric. 
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